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We report the first top quark mass measurement that uses a cross section constraint to improve the 
mass determination. This measurement is made with a dilepton tt event sample collected with the 

CDF II detector. From a data sample corresponding to an integrated luminosity of 1.2 fb^^ , wo mea- 
sure a top quark mass of 170.7l3;9(stat) ± 2.6(syst) ±2.4(theory) GeV/c^. The measurement with- 
out the cross section constraint results in a top quark mass of 169.7l| ;9(stat) ± 3.1(syst) GeV/c^. 

PACS numbers: 14.65.Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 
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The top quark is the heaviest known elementary par- 
ticle. Its mass Mt is a fundamental parameter in the 
standard model (SM). Together with the M^-boson mass, 
Mt places constraints on the SM Higgs boson mass 
At hadron colliders, the top quarks are mainly pair pro- 
duced via the strong interaction. Each top quark de- 
cays into a VF-boson and a fo-quark, and in the dilepton 
channel both T4/^-bosons decay to a charged lepton and a 
neutrino. The tt dilepton events have a small branching 
ratio, but they have a higher purity than single-lepton 
or all-hadronic final states. Because the two neutrinos in 
the final state are not detected, the dilepton channel top 
mass fit is under-constrained. However, measuring the 
mass in this channel is important because it provides an 
independent measurement of Mt that can be compared to 
measurements in other decay channels, allowing a consis- 
tency check of the ti hypothesis in the dilepton channel. 
Previous measurements of Mt in the dilepton channel are 
described in 2', 3, 4]. 

According to the SM, the theoretical tt cross section 
att has an exponential dependence on the top mass 0, @] . 
Therefore, the top mass value can be extracted from the 
observed event yield alone. By combining the theoreti- 
cal att{Mt) dependence with the top mass determination 
from the event kinematics, we can use the cross section 
information to improve the mass measurement, as re- 
ported in this Letter. 

In this novel measurement, the constraint provided by 
the mass dependent theoretical tt cross section is com- 
bined with a "template method" in which a top quark 
mass m^'^'^ is reconstructed for each event and in which 
the distribution of mj" '^'^ is compared with template distri- 
butions derived from simulation. In principle, any mea- 
surement of Mt and au could be combined with the the- 
oretical att{Mt) to obtain an improved top mass deter- 
mination. However, in our analysis we include the cross 
section constraint while properly taking into account the 
top mass dependence of the acceptance and all the cor- 
related systematic uncertainties. 

The template method adopted here is an enhanced ver- 
sion of the "full kinematic analysis" described in |2] . The 
enhanced version treats 6-tagged and non-tagged events 
separately. This separation improves the expected sta- 
tistical uncertainty by 20%; this represents a significant 
improvement over the previous analysis, which handled 
6-tagged and non-tagged events as a single sample. Intro- 
ducing the cross section constraint improves the expected 
statistical uncertainty further by 20%. In this paper, the 
measurement without the cross section constraint will be 
referred to as the "traditional" measurement. Both mea- 
surement techniques were fixed before the data were re- 
vealed. 

This measurement uses data collected by the CDF II 
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TABLE L Expected and observed number of events passing 
event selection criteria. Statistical and systematic uncertain- 
ties have been combined. 





Expected background 


diboson 


5.8 ± 0.9 


Z/7* ll,l = e, fi, T 


10.9 ± 2.3 


misidentified leptons 


8.8 ± 3.9 


Total 


25.6 ± 5.5 


tt {Mt = 170 GeV/c^) 


Expected signal 


62.1 ± 4.3 


Total expected 


87.7 ± 8.9 


Data 
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detector corresponding to an integrated luminosity of 
1.2 fb~^ The CDF II detector [tJ is a multi-purpose 
particle detector at the pp Tevatron Collider. Charged 
particle trajectories are measured with a silicon micro- 
strip detector and a drift chamber, which are immersed 
in a 1.4 T uniform magnetic field parallel to the beam 
directions. Electron, photon and hadron energies are 
measured with electromagnetic and hadronic calorime- 
ters. Muons are detected with drift chambers and scin- 
tillation counter hodoscopes located outside the calorime- 
ters. CDF employs cylindrical coordinates where 9 and 
(j) are the polar and azimuthal angles with respect to the 
proton beam. Transverse energy and momentum are de- 
fined as Et — Esm9 and px = psinO, where E is the 
energy and p is the momentum. 

The data for this analysis were collected using an in- 
clusive lepton trigger that required an electron or a muon 
with pt >18 GeV/c . After the offline reconstruction, 
the dilepton events were selected using the DIL selec- 
tion which requires two oppositely charged leptons 
with Pt >20 GeV/c Q, missing transverse energy 
due to the undetected neutrinos {^t > 25 GeV), and 
at least two jets with Et > 15 GeV. The main back- 
ground sources are diboson and Drell-Yan production 
and iy-|-jets events where a jet is mis-reconstructed as 
a lepton. Additional cuts were applied to reduce the 
background f^. The expected and observed numbers of 
events are summarized in Table HI After the event selec- 
tion, the sample was divided into two subsamples with 
significantly different signal-to-background ratios. The 
6-tagged sample includes 32 events in which at least one 
of the jets is identified as a 5-quark candidate through 
the presence of a displaced vertex [l^]. This subsample 
has an expected signal-to-background ratio of 11:1. The 
non-tagged sample comprises 45 events in which none of 
the jets is identified as a 6-quark candidate. In this sub- 
sample the expected signal-to-background ratio is 1:1. 

Because the two neutrinos are not detected, the recon- 
struction of the top quark mass from dilepton events is 
under-constrained. Top mass reconstruction can be ac- 
complished by considering a kinematic variable which is 
not observable on an event-by-event basis, but the distri- 



5 



bution of which is predictable and independent of the top 
mass value. In this analysis the distribution of p", the 
longitudinal momentum of the tt system, was adopted as 
the variable. Monte Carlo simulations, generated with 
PYTHIA and the CDF II detector simulation [iJI , in- 
dicate that the distribution of is nearly independent 
of the top mass, and is described by a Gaussian distribu- 
tion with a mean of zero and width of 195 GeV/ c. The 
validity of the Monte Carlo simulation was tested with 
data from the lepton-|-jets decay channel where p** can 
be explicitly reconstructed. 

For each event, a top mass mj''^^ is reconstructed from 
the event kinematics as follows. The jet energies are cor- 
rected for defects in the calorimeter response, computed 
from Monte Carlo simulations, to correspond to the en- 
ergies of the primary 5-quarks. After these jet energy 
corrections, the two components of transverse missing 
momentum are taken as the sum of the neutrino trans- 
verse momentum components. Along with assumptions 
on the masses of the final state particles and additional 
constraints on Mw± = 80.4 GeV/c^, Mt = Mi, and 
Pz +P* = Pz*! ^ top mass can be calculated 0]. 

A wide range of possible values is incorporated by 
calculating the top mass 10000 times. For each itera- 
tion, p** is randomly drawn from its expected distribu- 
tion. Similarly, the jet energies and are smeared ac- 
cording to their resolutions. For each iteration, if a so- 
lution is not found using the fixed values of My/ and 
Mt, solutions within Mw± = 80.4 ± 3.0 GeV/c^ and 
Mt = Mt± 2.0 GeV/c^ are accepted. 

For a given event, we obtain two distributions of pos- 
sible top quark masses, each corresponding to a different 
lepton-jet pairing. The pairing which has the higher frac- 
tion of solutions is selected. This choice is correct for 70% 
of simulated tt events. If the number of entries in this 
distribution is less than 100, the event is rejected. Ac- 
cording to Monte Carlo studies, 91% of signal and 78% 
of background events pass this mass reconstruction re- 
quirement. The most probable value of a spline fit to 
the distribution selected is taken as a per-event top mass 
ml"". 

Templates of reconstructed top mass distributions were 
created from various ti and background samples. Sig- 
nal templates were generated from ti Monte Carlo sam- 
ples with generated top masses ranging from 150 to 
200 GeV/c^, separately for 6-tagged and non-tagged sig- 
nal events. Diboson and Z —>■ II templates were gen- 
erated from Monte Carlo simulation. A template for 
misidentified leptons was created using data. The back- 
ground templates were combined according to the ex- 
pected contribution of each background source. It was 
observed from simulation that using the same common 
background template for 6-tagged and non-tagged sam- 
ples provides as good a performance as using separate 
templates. The common background template was there- 
fore used for both subsamples. 

In the traditional measurement, the top mass is ex- 
tracted by comparing the reconstructed mass distribu- 



tions from data to the signal and background template 
parametrizations (/« and /{,, respectively) using an un- 
binned likelihood fit. The likelihood includes free param- 
eters for the number of signal events rig and background 
events nb in each subsample, and for the top mass Mt- 
The total likelihood takes the form 



where each of the subsample likelihoods is: 



(1) 



-'non/6— tagged 
■^shapc 

r 

- In Chg 



TV rig x/s(m[^'',Mt)+nb x/b(m^^'') 



n 



e''"'' + "b'(«,+nb)^ 

("b-"b ) 

■ 



(2) 



In the formula, TV is the number of data events in the sub- 
sample, n,^^ and are the expected number of back- 
ground events and its uncertainty, and m^'r" is the value 
of mil"" for the i'th event. The top quark mass hypothesis 
which minimizes -ln(£) is taken. 

To test the method, we performed Monte Carlo exper- 
iments of signal and background events. The numbers 
of signal and background events in each experiment were 
varied according to Poisson distributions using the ex- 
pected mean numbers of events. According to the Monte 
Carlo experiments, the method is unbiased and returns 
appropriate uncertainties. 

In 1.2 fb^^ of data, 31 6-tagged and 39 non-tagged 
events pass the event selection criteria and have so- 
lutions for m\"" . The mass measurement from the 
6-tagged sample gives Mt = 171.9 j:g:4(stat.) GeV/c^ 
and the measurement from the non-tagged sample gives 
Mt = 166.0t7 g(stat.) GeV/c^. Applying the traditional 
method to the two subsamples with a joint likelihood, 
we measure Mt = 169.7j:|:9(stat.) GeV/c^. In 34% of 
Monte Carlo experiments, the statistical uncertainty was 
smaller than the measured statistical uncertainty from 
data. The reconstructed top mass distribution from data 
is shown in Fig. [1] 

The top mass measurement can be further improved 
by taking into account the top mass dependence of the tt 
production cross section. The expected number of signal 
events can be expressed as 



ns{Mt) ^ ati{Mt) ■ a{Mt) ■ L ■ pl 



(3) 



where ati{Mt) is the theoretical tt cross section, a{Mt) is 
the acceptance of ti events, L is the integrated luminosity, 
and Pmass is the probability of obtaining a solution for 
ml"". 

The dominant top mass dependence on the expected 
number of events comes from the theoretical att- We use 
a NLO calculation of (Jtt evaluated at three different top 
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Reconstructed Mass ( GeV/c^ ) 

FIG. 1: Reconstructed top mass distribution from data to- 
gether with the signal and background parametrizations. 



masses [5|; we parametrize the mass dependence of (Tjt 
on the top mass using the functional form described in 

atiiMt) - 6.70 • e("5-M,)/32.29 (4) 

The acceptance a{Mt) was studied using tt Monte Carlo 
simulation, separately for 6-tagged and non-tagged sam- 
ples. The Monte Carlo acceptances were corrected for 
trigger efficiencies and for scale factors arising from dif- 
ferences between data and simulation. The combined 
Monte Carlo acceptance corrections are between 74% and 
95%, depending on the lepton flavor and pseudorapid- 
ity. The dependence of the acceptance on the top mass 
is linear, increasing about 30% in the top mass range 
of 150 to 200 GeV/c^. The integrated luminosity, L, is 
1118 pb^^ for the 6-tagged sample and 1189 pb^^ for the 
non-tagged sample. The two integrated luminosities are 
different because of the smaller data set using the sili- 
con detector, which is required for 6-tagging. The signal 
mass reconstruction probability, p^ass ^ "^^^ measured to 
be 91±1.1% for both 6-tagged and non-tagged samples, 
and was found to be independent of the top mass. 

The cross section constrained top mass measurement 
uses information from the reconstructed top mass distri- 
bution as well as the observed number of events. The 
per-event mass reconstruction method and the template 
parametrizations are exactly the same as in the tradi- 
tional measurement. The information from the num- 
ber of events is added in the likelihood function by re- 
placing Us in Eq. ([2]) with ns{Mt) from Eq. ([3]); thus 

C = /:b-taggcd(Mt,?l[^) X /:„on-taggcd(Mt,nr")- The 

number of background events nb and the top mass Mt 
are free fit parameters as in the likelihood function of the 
traditional measurement. The uncertainty in the theoret- 
ical modeling of Uti is not included in the likelihood; it is 
treated in the same way as other systematic uncertainties 
described below. 



Simulated experiments are used to verify that the 
cross section constrained method is unbiased and re- 
turns appropriate uncertainties. We measure Mt — 
170.7!t3 g(stat.) GeV/c^. The statistical uncertainty is 
consistent with expectations: 34% of Monte Carlo ex- 
periments returned an uncertainty smaller than the un- 
certainty measured from data. 

The sources of systematic uncertainties are summa- 
rized in Table |TT1 The jet energy scale uncertainty is 
dominated by the uncertainty in jet energy corrections. 
This uncertainty was studied by shifting the jet energies 
by ±lcr, and half of the mass difference was taken as the 
systematic uncertainty. Since the jet energy corrections 
were determined for light quark jets, we evaluated an ad- 
ditional systematic uncertainty from possible differences 
between b jets and light quark jets [ij]. The total un- 
certainty from the jet energy scale is 1.8 (2.9) GeV/c^ 
for the cross section constrained (traditional) measure- 
ment. The cross section constrained measurement is less 
sensitive to the jet energy corrections because a change 
in the jet energy scale shifts the top mass determina- 
tion from the event yield in the opposite direction to 
that from kinematic reconstruction. The signal modeling 
uncertainty is 0.9 (0.8) GeV/c^, and takes into account 
differences in parton showering between the PYTHIA [ii| 
and HERWIG Monte Carlo generators, uncertainties in 
initial and final state radiation modeling, and differences 
in parton distribution functions between MRST [l^ and 
the full set of CTEQ6M [l^ eigenvectors. Possible im- 
perfections in modeling the Z ^ II and misidentified 
lepton backgrounds combine to give 0.3 (0.3) GeV/c^ 
background modeling uncertainty. The uncertainty from 
template statistics is 0.4 (0.5) GeV/c^. A 1% uncer- 
tainty in the lepton px introduces an uncertainty of 
0.2 (0.2) GeV/c^. The cross section constrained measure- 
ment has an additional uncertainty of 1.6 GeV/c^ from 
the expected number of events. This uncertainty includes 
1.1 GeV/c^ uncertainty from the integrated luminosity, 
0.5 GeV/c^ from the acceptance, 0.9 GeV/c^ from the 
expected number of background events and 0.5 GeV/c^ 
from the mass reconstruction probability. 

TABLE II: Summary of systematic uncertainties for the tra- 
ditional (T) and cross section constrained (C) measurements. 



Systematic Source 


AMt (GeV/c") 




T 


C 


Jet energy scale 


2.9 


1.8 


Signal modeling 


0.8 


0.9 


Background modeling 


0.3 


0.3 


Template statistics 


0.5 


0.4 


Lepton pt 


0.2 


0.2 


Expected number of events 


n.a 


1.6 


Total 


3.1 


2.6 



The uncertainty in the theoretical att{Mt = 
175 GeV/c2) is toll pb @. We propagated this un- 
certainty to the top mass by changing the number of 
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signal events in the Monte Carlo experiments. The es- 
timated uncertainty on the top mass is 2.4 GeV/c^. 
Simulation studies show that this cross section con- 
strained top mass measurement is not very sensitive to 
the probability shape of the theoretical att uncertainty. 
Figure [2] shows the cross section constrained top mass 
measurement in the Mt - Uti plane. The extracted 
top mass from the cross section measurement p7j only 
is 178.3l8°o^(exp)t^;°(theory) GeV/c^, consistent within 
about one standard deviation with the result from the 
traditional analysis. 
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FIG. 2: The measured cross section constrained top mass 
is shown in the Mt - att plane. The innermost error bars 
correspond to the statistical uncertainty, the middle ones 
the statistical-|-systematic uncertainty, and the outermost er- 
ror bars show the statistical-|-systematic-|-theory uncertainty. 
The hatched areas mark the traditional top mass measure- 
ment and the measurement in the dilepton channel with 
statistical-|-systematic uncertainties. 



In summary, we have introduced a new way to improve 
the template top mass measurement in the dilepton chan- 
nel by using a theoretical cross section constraint. With 
this measurement, we compare the reconstructed top 
mass distribution to templates and the observed number 
of events to expectation. In 1.2 fb^^ of data collected by 
the CDF II detector, we measure a top quark mass of 
170.7t3:g(stat) ± 2.6(syst) ± 2.4(theory) GeV/c^. This 
measurement is in good agreement with the top 
mass measurement made without a cross section con- 
straint, which gives 169.7t|:9(stat) ± 3.1(syst) GeV/c^ 
and with top quark mass measurements made in other 
decay channels [H, [H, [l^l . 
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